How to detect gap nodes of a superconductor by angle-resolved specific heat 

experiment? 
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The specific heat oscillation in the mixed state of type II superconductors is studied theoretically 
when rotating field within a plane containing gap minimum and maximum. The calculations are 
performed microscopically by solving quasi-classical Eilenberger equation for vortex lattices. The 
field dependence of the oscillation amplitude can discriminate between the nodal and anisotropic 
gap with a finite minimum and the oscillation phase gives the gap minimum position on the Fermi 
surface. These also provide a way to separate out the anisotropic behavior due to the Fermi velocity. 

PACS numbers: 74.25Bt, 74.250p, 74.20Rp 



I. INTRODUCTION 

There has been much attention focused on exotic su- 
perconductors, including high Tc cuprates and heavy 
Fermion materials in recent years. In addition to the 
spin structure or parity of the Cooper pair, the orbital 
function or the gap structure on the Fermi surface is deci- 
sive to characterize its superconductivity. These studies 
are expected to lead to a new pairing mechanism. Even 
in conventional superconductors the energy gap can vary, 
depending on the position on the Fermi surface. The de- 
gree of the anisotropy in the gap function is an important 
factor in understanding a superconductor in question. To 
distinguish a nodal superconductor from an anisotropic 
one with a finite minimum gap is of particular impor- 
tance because the gap function can not change its sign in 
the latter while it can in the former. Also it is crucial to 
determine the maximum and minimum gap positions on 
the Fermi surface. 

It is now widely recognized that the zero-energy den- 
sity of states (ZEDOS) sensitively reflects the gap struc- 
ture, which is probed by a variety of experimental meth- 
ods such as specific heat, thermal conductivity or scan- 
ning tunneling spectroscopy. This is particularly true for 
physical quantities in the mixed state of type II supercon- 
ductors. Induced vortices under an applied field carry a 
certain amount of ZEDOS around each vortex core which 
depends on the gap structur e^ i^i'^ . The Sommerfeld coef- 
ficient 7(i?) as a function of magnetic induction B, which 
is nothing but ZEDOS induced by vortices, is found to be 
one of physical quantities to reveal the gap topology. In 
fact, there have been several ^{B) experimenta^i^iSiLSiS 
on such as 2H-NbSe2, VsSi, NbsSn or CeRu2. We re- 
cently demonstrate that precise measurement ^{B) at 
low field gives rise to indispensable information on the 
gap anisotropy In order to better characterize the gap 



structure it is urgent to provide further ways to analyze 
experimental data. For example, anisotropic behaviors in 
a superconductor could come from the two main sources. 
One is the gap structure itself and the other is the Fermi 
velocity anisotropy due to band structure. It is often the 
case that these two kinds of anisotropy are mixed up and 
difficult to separate out individually, leading to an am- 
biguous conclusion as for the gap structure. Thus we are 
required to devise some method to disentangle these two 
anisotropy effects. 

Recently the angle-dependent specific heat experi- 
ments for the mixed state in several superconductors 
LuNiaBsC El CeCoIng and Sr2Ru04 ,13 have been 
performed to yield characteristic oscillation pattern in 
^{B). A few percent oscillation amplitude relative to the 
total in these experiments is generally consistent with 
the theoretical estimateS. for nodal superconductors or 
strongly anisotropic gap superconductors. However, it 
remains open to distinguish between them. Specifically, 
it is reported that the oscillation amplitude becomes van- 
ishing toward _B — s- in Sr2Ru04 .13. while it remains a 
finite value in LuNi2B2C [TJ and CeCoIns 0[ It was 
speculated that in the former either the gap structure 
has a finite minimum gap or the gap in the minor band 
may mask the oscillation in lower fields. Thus we need 
to know precise behavior of the field-dependence of the 
oscillation amplitude for the two cases. 

The purpose of the present paper is to examine the os- 
cillation amplitude of the angle-dependent Sommerfeld 
coefficient ^{B) when B rotates within a plane contain- 
ing the gap minimum and maximum for several typical 
gap topologies, including line and point node supercon- 
ductors and a superconductor with a finite minimum 
gap. We also study the oscillation behavior of "i{B) for 
isotropic gap case with the anisotropic Fermi velocity. It 
turns out that this anisotropy also yields a substantial 
specific heat oscillation under field rotation, but we will 
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provide information on how to distinguish it from the gap 
anisotropy case. The existing data of the j{B) oscilla- 
tion experiments on LuNi2B2C, CeCoIns and Sr2Ru04 
are analyzed in this regard. After a brief introduction to 
quasi-classical framework to show how to calculate the 
ZEDOS for various situations, we describe the results of 
the ZEDOS oscillations to examine the differences be- 
tween node versus nodeless gap cases in Section 3 and 
also show the results for anisotropic Fermi velocity in 
Section 4. The Section 5 is devoted to the point node 
case in comparison with the line node case to supplement 
the above analysis. In the Section 6 we give a summary 
and discussions on the existing data. All computations 
are done assuming a Fermi sphere. 



II. QUASI-CLASSICAL THEORY AND ZEDOS 



Here, Jpg is integral over the Fermi surface, lud is cut off 
frequency, Nq is the total density of states for one spin 
at the Fermi surface in the normal state 



Nn = 



(5) 



FS 



and p{kp) is the angle resolved density of states at the 
Fermi surface: 



pikp) = 



1 1 



(27r)37Vo \vy 



normalized so that: 



J (fkFpikp) = 1. 

FS 



(6) 



(7) 



The amplitude of the specific heat oscillations ap- Density of states N{r,E) is defined as 



.11.12.13 



as men- 



pears to be very small, just few percent; 
tioned before. The smallness of the effect necessitates 
the use of numerical solutions of Gorkov's microscopic 
equations of superconductivity to accurately estimate the 
specific heat amplitude. Quasi-classical approximation 
of Gorkov's equations, that we solve numerically here, 
is good as long as condition kpS, » 1 is met. Here, 
kp^ is in the order of atomic length, and ^ is coherence 
length. Pairing interaction between electrons is modeled 
as V{v,v') = Voil,{v)rt{v'). This greatly simplifies the 
analysis, since pairing potential in this model can be 
written as A{r,v) ~ n{v)^>{r). The orbital part of the 
pairing potential, r2(v), we simply call: the gap function. 
Eilenberger equations then read as (fi = 1) 



[2uj + vU] f{cu, r, v) = 2A(r, v)g{cu, r, v), (1) 



[2co - vU*] f{co, r, v) = 2A*(r, v)g{u, r, v). (2) 

Here 11 = V + (27ri/$o)^ is gauge invariant gradient, 
A is vector potential, $0 is flux quantum, v — v{kp) 
is Fermi velocity defined as v(k) = WkE{k), with E{k) 
being energy function of electrons in band; Fermi wave 
vector kp can be found from the equation Eik) — Ep; 
Lo = 7rT(2n -I- 1) with integer n is Matsubara frequency, 
g and / are normal and anomalous Green's function and 
p{Lo,r,v) = f*{io,r,—v). Normalization condition for 
the Green's functions is g"^ + ff^ = 1. 

The selfconsistency equations for the gap function and 
current density are: 



A(r, v) = 2^iVoT ^ / d^k'p Viv, v')p{k'p)f, 
j = 47ri|e|iVor^ / (fkp p{kp)vg. 

CJ>0 pg 



(3) 



(4) 



N{r,E) ^ No J SkpKe g{iuj ^ E + i6,r,v)p{kp) (8) 

FS 

We are interested in low temperature zero-energy den- 
sity of states (ZEDOS) N{r,E = 0). This is because in 
the limit of small temperatures, T — > 0, ratio of specific 
heat in superconducting state Cg and normal state C„ is 
given by 



^.^C^^Nir,E = 0) 

T^O Cn No 



(9) 



Here N{r, E = 0) is spatially averaged ZEDOS in su- 
perconducting state. In our calculation we set T ~ 
Tc/10. Since ZEDOS and specific heat are proportional 
at low temperatures we use these two terms concur- 
rently throughout the text. Numerical procedure for 
solving quasi-classical equations of superconductivity is 
described in Ref. Magnetic field is measured in units 
<i>o/(27ri?o), where $0 is flux quantum, and Rq — 0.882 ,^0 
(^0 is BCS coherence length). 

Angular dependence of ZEDOS is already studied nu- 
merically for some typical cases of nodal gap function: 
3D d-wave, polar state, and axial stated. Also, ZEDOS 
is studied by using the high field approximate solution 
of Eilenberger equationsi^tiSiiSiiii^ The focus in that 
study was on field dependence of ZEDOS for characteris- 
tic magnetic field directions. Just to recall, as one rotates 
magnetic field from the gap node direction toward gap 
maximum direction, ZEDOS increases in low field, while 
it decreases in fields near i?c2- lu other words, specific 
heat oscillation amplitude changes sign with increasing 
field. In this paper we compare ZEDOS oscillation am- 
plitude for nodal and nodeless superconductors in the 
limit of low fields. Along with the nodal gap structures 
already studied in Ref. 3, we present also results for some 
other typical cases known in the literature. 
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III. NODAL GAP VS NODELESS GAP 

In order to distinguish the nodal gap superconductor 
from the nodeless superconductor with a finite minimum 
gap, we examine the following model for the gap structure 



il{(p,9) — r^o (a) V 1 + a cos iip 



(10) 



with if polar angle and 9 azimuthal angle in polar co- 
ordinates. The parameter a measures the degree of 
anisotropy. For a = the gap function is isotropic, while 
for a = 1 the gap function reduces to two-dimensional 
(2D) version of d^^-y^ -wave superconductivity. We 
choose the prefactor 17o(a) so that the average of fi^ over 
the Fermi surface is unity, independent of a. ZEDOS 
is calculated for fields rotating in the basal plane with 
= 7r/2. Only two field directions are of our interest, B 
along the gap minimum ((^ (2fc + \)tt/A), and B along 
the gap maximum (</? = /c7r/2). Here, k is integer. We 
performed the calculation for two anisotropy parameters: 
a = 0.5 and a = 1. In both cases, low field ZEDOS is al- 
ways minimum when the magnetic field is oriented along 
the gap minimum (node) and vice versa. 

In Fig. Hthe ratio R = N{E = 0, = 0)/N{E = 
0,ip = j) of maximum and minimum ZEDOS, for field 
rotating in the basal plane, is plotted for both anisotropy 
parameters. As is clear from Fig. ^ there is a striking dif- 
ference in low field dependence of ratio R for nodal a = 1 
and nodeless a = 0.5 superconducting gap. Namely, for 
the nodal gap case (a = 1) there is a finite amplitude 
ZEDOS osciUation {R ^ 1) towards B ^ 0. This is 
contrasted with the nodeless gap case (a — 0.5) where 
R becomes unity at lower fields, showing a maximum 
at the intermediate field region B„iax- This means that 
the specific heat oscillation diminishes there. The field 



Braax comcs from the physical reason that at i? < Bmax 
the spatial extension of the ZEDOS is confined to each 
vortex core region, yielding more or less isotropic ZE- 
DOS landscape. This ZEDOS feature does not cause the 
specific heat oscillation in this lower field. The spatial 
extension of the ZEDOS landscape depends on the size 
of the minimum gap because the gap acts as a potential 
for quasi-particles, that is, the zero-energy quasi-particles 
are strongly confined and localized near each core. Thus 
Bmax signals the characteristic field where the localized 
zero-energy quasi-particles begin overlapping and tends 
to become smaller as the minimum gap decreases. Since 
in the nodal gap case the spatial extension of the ZEDOS 
is extended to infinity, Bmax approaches zero, indicating 
that R stays a constant towards smaller fields. 



1.02 



1.01- 
1.00- 
0.99- 
0.98- 
CC 0.97- 
0.96- 
0.95- 
0.94- 
0.93- 



j\nisotropic Fermi velocity 



0.00 



0.20 



0.40 0.60 
B 



0.80 



1.00 



FIG. 2; Ratio R = N{E = 0, ^ = 0)/N{E = 0,^ = f) of 
minimum and maximum ZEDOS in case of anisotropic Fermi 
velocity modeled by Eq. IIIH . 
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FIG. 1: Ratio 7? = N{E = 0,(p = 0)/N{E = 0, <^ = f) 
of minimum and maximum ZEDOS for fields rotating in the 
basal plane of the crystal. Full circles are for anisotropy pa- 
rameter a = 0.5, while full squares are for 2D d-wave super- 
conductor (o = 1). 



IV. ANISOTROPIC FERMI VELOCITY 

We analyze the other type of anisotropy. The super- 
conducting gap is assumed to be isotropic, or to have the 
same value all over the Fermi surface, but the Fermi ve- 
locity itself is anisotropic. The amplitude of the Fermi 
velocity on the Fermi sphere is modeled as: 

v{ip,e) ^vo{b){l + bcos4ip). (11) 

The parameter b measures the degree of anisotropy. For 
convenience the prefactor vo{b) is chosen so that den- 
sity of states in normal state is independent of b. Four- 
fold variation of the Fermi velocity in the basal plane 
is a quite simple model but will suffice for our purpose. 
Even in this hypothetical case with the isotropic gap and 
the anisotropic Fermi velocity, ZEDOS, and thus specific 
heat, depend on magnetic field direction. In Fig. 13 we 
plot the ratio R of minimum and maximum ZEDOS as 
a function of magnetic field. The anisotropy parameter 
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b = 0.5 is used in the calculation. The ratio R between 
maximum and minimum ZEDOS is of the order of a few 
percent. This is the same order of magnitude as in the 
above gap anisotropy cases. This is also the same or- 
der of magnitude observed^^in LuNi2B2C. This warns us 
that the effect of the Fermi surface structure on the spe- 
cific heat oscillation cannot be neglected in intermediate 
fields. However, it is seen from Fig. 2 that extrapo- 
lation of maximum/minimum ratio R to B — gives 
R = I, i.e. disappearance of specific heat oscillation. 
This is qualitatively the same behavior as seen in the 
anisotropic gap case with a finite minimum gap, shown 
in Fig. 1. Physically this is due to the spatial extension 
of the zero-energy quasi-particle state: In the Fermi ve- 
locity anisotropy case ZEDOS is confined to each core in 
a rather isotropic gap manner. Thus to induce the oscil- 
lation finite field is needed, above which because of the 
overlapping of the zero-energy quasi-particles, ZEDOS 
begins to exhibit an oscillation. 

It is also noticed that in the above 6 > case R < 1, 
meaning that the four-fold oscillation pattern is maxi- 
mally phase-shifted by 7r/4 from the gap anisotropy cases 
with i? > 1 in Fig. ^ Needless to say, the sign b is arbi- 
trary for a given material, but it is physically plausible 
case 6 > when a > because in the angle-resolved DOS 
N((p) cx l/v{ip) the larger energy gap {(p = 0) coincides 
with the larger angle-resolved DOS. This is indeed the 
case in boro-carbides. Therefore, we can clearly distin- 
guish the two anisotropy effects by measuring the angle- 
dependent specific heat to monitor both the oscillation 
amplitude and its phase. 



V. POINT NODE GAP 

In this section we consider the point node gap struc- 
ture. Among several possible point node topologies we 
take up typical examples, the so-called axial state and 
"s-|-g" model. The former is known to be realized in su- 
perfluid "^Hc A phase and the latter is an candidate for 
boro-carbides22,. We also consider the polar state with a 
line node for comparison. 



A. Axial and polar gap function 

The polar gap function has a horizontal line node in 
the crystal basal plane, while the axial gap function is 
characterized by two point-like nodes at the poles of the 
Fermi sphere. In polar coordinates with (p and 6 denote 
polar and azimuthal angle, then the polar gap function 
is presented with 

n{ip,e) = Vscose, (12) 

and the axial gap function is presented with 

n{(p,9) ^ y^sinO. (13) 
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FIG. 3: Ratio R = N{E = O,antinodal)/N{E = 0, nodal) as 
a function magnetic induction for polar and axial gap func- 
tion. 
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FIG. 4; a) Field dependence of ZEDOS for two field direc- 
tions: nodal and antinodal when fields are rotated in the basal 
plane, b) Ratio R = N{E = 0, antinodal) /N{E = 0, nodal) 
as a function magnetic induction. 



The low temperature ZEDOS for these two gap func- 
tions is studied in detail in Ref. ^ By rotating mag- 
netic field in a plane that contains c-axis, ZEDOS peri- 
odically changes. The ZEDOS maximum appears for the 
field oriented along the gap node, and the ZEDOS min- 
imum appears when magnetic field is directed along the 
gap maximum. Although field dependence of ZEDOS 
is already presented in Ref. 0, for the purpose of this 
paper we present those data in a slightly different form. 
Namely, we are interested in the amplitude of ZEDOS os- 
cillation as a function of B. Therefore, in Fig. |3|the field 
dependence of ratio R — N{E = , antinodal) / N [E = 
0, nodal) is shown. The low field ratio differs significantly 
for the axial and polar gap functions. The important 
point to notice is that in both cases ratio i?, in the limit 
of low fields, differs from 1. This means that there is a 
finite amplitude of ZEDOS oscillation as H — > 0. This 
is in accord with the conclusion in the previous sections. 
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B. "s+g" model 

It has been argued that the gap function in boro- 
carbides is a mixture of s-wave and g-wsve supercon- 
ductivity with aptly chosen weighting factors so that the 
pairing potential 9) does not change sign on the 
Fermi surface but has point-like zeros: 

f^(^,^) = ^^(l-sin4 0cos4^). (14) 

Like in previous examples of anisotropic pairing func- 
tions, we choose the Fermi surface as a sphere to see the 
effect of the gap structure on low temperature thermody- 
namics. As we have already shown, even the anisotropy 
of Fermi surface alone can account for direction depen- 
dent specific heat. The effect of Fermi surface disappears 
only in the limit of low fields. Having in mind that band 
structure in boro-carbides is far from being isotropic, this 
simple "s-|-g" gap function on Fermi sphere, may not be 
an appropriate model which can accurately estimate the 
amplitude of specific heat oscillation in boro-carbides, 
but we can gain qualitative tendency in this case. In 
Fig. ^) the field dependence of low temperature zero- 
energy DOS is shown for nodal and antinodal field direc- 
tions. We define the antinodal direction as a direction in 
the basal plane with maximum value of superconducting 
gap. Magnetic induction B is scaled with i7c2 (antinodal) 
and Hc2 (nodal) respectively. Note that Hc2 (antinodal) > 
ifc2 (nodal). As one may already anticipate, difference 
between the maximum ZEDOS (field along antinodal di- 
rection) and the minimum ZEDOS (field along nodal di- 
rection) remains finite in low fields. Thus their ratio 
N{E = 0, antinodal) /N{E = 0, nodal) ^ 1. This is 
shown in Fig^). Noteworthy is comparison with gap 
structure of the axial state fl{ip,9) — ■\/3/2sin6'. Both, 
the axial state and "s-l-g" model, have point like nodes, 
while their field dependence of ZEDOS is different. The 
7(-B) behavior as a function of B is not governed only by 
the gap node topology, point-like or line- like nodes. It is 
not unique for all types of gap structures with point-like 
(or line-like) nodes. It rather reflects the gap value on 
average, and it is rather sensitive to the functional form 
of the gap function in the vicinity of node. 

VI. DISCUSSIONS ON THE SPECIFIC HEAT 
EXPERIMENTS 

In this section we discuss the angle dependent spe- 
cific heat experiments on three materials, LuNi2B2C, 
CeCoIns and Sr2Ru04 in the light of the present cal- 
culations. 

A. Boro-carbides 

Park et afcitepark measure the angle-dependent j{B) 
for LuNi2B2C and detect the four-fold oscillations in var- 



ious fields, identifying the gap minimum is located in 
[100] direction of tetragonal crystal because the oscilla- 
tion maximum is in [110] direction. This identification 
is supported by our present result. The oscillation am- 
plitude becomes smaller as B decreases. However, it is 
rather difficult to judge whether the gap has a node or 
a finite minimum gap from their experiment where the 
detailed low field data are lacking. In connection with 
other experiments which suggest strong anisotropic gapSi 
or point node gap22, in this system, it is interesting to ex- 
tend their measurement to lower field to determine the 
precise gap structure. 

Izawa et aS^ show in their angle-dependent thermal 
conductivity measurement: (1) In the similar oscillation 
pattern the maximum appears in [110], coinciding with 
that in Park et al. (2) The oscillation amplitude dimin- 
ishes when field rotates conically out of the basal plane. 
For the polar angle 6 > 7r/4 it almost vanishes. They 
conclude a point node gap located along [100] direction. 
This assertion is based on a theoretical calculationSS of 
angle dependence of thermal conductivity. It may be in- 
formative in this connection to show our result: We have 
also performed the calculationp2^ when field rotates con- 
ically to check how the oscillation amplitude varies as 
a function of the polar angle 9. It is found that it de- 
creases quickly as 9 increases from zero for both vertical 
line node and the "s-|-g" point node cases, being unable 
to distinguish these two cases by conical field rotation of 
specific heat experiment. 



B. CeCoIng 

Aoki et ai^ perform the angle dependent specific heat 
experiment in this system, observing a substantial os- 
cillation amplitude when rotating field within the basal 
plane of tetragonal symmetry crystal. In the oscillation 
pattern the maximum occurs for [110] direction. This 
suggests the d^y gap function because the oscillation am- 
plitude stays constant towards lowest fields. This conclu- 
sion appears to be inconsistent with the angle-dependent 
thermal conductivity experiment by Izawa et aS^, who 
conclude the dx2-.y2 gap function. Note, however, that 
their data themselves exhibit the oscillation maximum 
for [110] direction, consistent with Aoki et a&^. 

C. Sr2Ru04 

According to Deguchi et al^^, who measure the angle- 
dependent specific heat on this system by rotating field 
within the basal plane of tetragonal crystal, the four-fold 
oscillation amplitude decreases below a threshold field 
~ 0.3T and changes its sign near Bc2 ^ 1.5T. The ex- 
istence of the threshold field is in accord with our cal- 
culation where the gap structure has a finite minimum 
gap, definitely excluding the vertical line node in the so- 
called main 7 band. (We can not say anything about 
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the possible horizontal line node^.) An interesting point 
in this system is the fact that we know accurately the 
Fermi velocity anisotropy in the 7 band, where the larger 
Fermi velocity is directed to [110]. Since the observed 
oscillation maximum occurs along [110], the experiment 
unambiguously excludes the oscillation due to the Fermi 
velocity anisotropy. Combining these two facts we con- 
clude that there exists the anisotropic gap structure with 
a finite minimum in the basal plane of the main 7 band in 
Sr2Ru04. We can not commit ourselves on further con- 
clusion concerning the minor a and /3 bands or the hor- 
izontal line node based on the existing data by Deguchi 
et aH^. 



VII. SUMMARY AND CONCLUSION 

In this paper, we have calculated the zero-energy den- 
sity of states in the mixed state at low temperature by 
employing quasi-classical Eilenberger formalism, which 
is valid for a wide variety of superconductors. We have 
focused on the angle dependence of the zero-energy den- 
sity of states, which is directly measured through specific 
heat experiment as the Sommerfeld coefficient, for super- 
conductors with both nodal or nodeless gap structures. 

We have demonstrated that the specific heat angular 
dependence provides useful information concerning the 
gap structure, namely, the position of the node or the 



gap minimum on the Fermi surface and also the existence 
or non-existence of the gap node. Furthermore, we give 
information to distinguish two sources of the anisotropy, 
either due to the gap itself or due to the Fermi velocity 
of band structure. These proposed methods, we believe, 
add yet another dimension to firmly establish the gap 
structure. 

A few examples studied here does not exhaust all pos- 
sibilities for superconducting gap function. Neither all 
kind of different Fermi surface structures can be covered. 
Nevertheless, from the comparison of field dependence 
of ZEDOS in: a) nodal, and b) fully gaped anisotropic 
superconductors, one can conjecture a behavior that is 
common for each group of superconductors. In nodal 
superconductors the specific heat oscillation amplitude 
persists down to very low fields. In contrast to this be- 
havior, fully gapped superconductors, as one decreases 
magnetic field the oscillation amplitude gradually dimin- 
ishes. To remind the reader again, only low temperature 
specific heat is discussed here. 
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